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Abstract. We present a study of the stellar parameters, distances and spectral energy distributions (SEDs) of 
HD 34282 and HD 141569, two pre-main sequence Herbig AeBe stars. Both objects have been reported to show 
'anomalous positions' in the HR diagram in the sense that they appear below the main sequence. A significant 
result of this work is that both stars are metal-deficient. The Hipparcos distance of HD 34282 is very uncertain 
and the current study places the star at the expected evolutionary position in the HR diagram, i.e. as a PMS star. 
The distance for HD 141569 found in this work matches the Hipparcos distance, and the problem of its anomalous 
position is solved as a result of the low metallicity of the object: using the right metallicity tracks, the star is in 
the PMS region. The SEDs are constructed using data covering ultraviolet to millimetre wavelengths. Physical, 
non-parametric models, have been applied in order to extract some properties of the disks surrounding the stars. 
The disk around HD 34282 is accreting actively, it is massive and presents large grains in the mid-plane and small 
grains in the surface. HD 141569 has a very low mass disk, which is in an intermediate stage towards a debris-type 
disk. 
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1. Introduction 

The study of protoplanetary disks is currently undergo- 
ing an exciting stage partly propelled by the discovery of 
extrasolar planetary systems following the detection of 51 
Peg B by Mayor & Queloz (1995). The possible discovery 
of telluric planets in the near future, in addition to the 
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jovian-like planets already detected, will pose interesting 
questions on the formation of extrasolar planetary sys- 
tems. Knowledge of the properties of protoplanetary disks 
and how they evolve to debris disks around main-sequence 
(MS) stars is one of the tools required for modelling such 
a process. 

Observed spectral energy distributions (SEDs) of pre- 
main sequence (PMS) stars are widely used to study the 
properties of protoplanetary disks and to classify T Tauri 
and Herbig AeBe stars (HAeBe hereafter) into an evo- 
lutionary scheme (Adams et al. 1987; Hillenbrand et al. 
1992). Although the interpretation of a given SED is, to 
some extent, model dependent, the theoretical modelling 
of SEDs constitutes an invaluable tool for understanding 
the structure and properties of protoplanetary disks, e.g. 
Chiang & Goldreich (1997, 1999), D'Alessio et al. (1998, 
1999, 2001), Dullemond et al. (2001). The EXPORT con- 
sortium (Eiroa et al. 2000) observed a large sample of PMS 
and Vega-type stars during the 1998 International Time 
Programme of the Canary Islands' Observatories. One of 
the driving goals of this effort was the study of the evolu- 
tion and properties of protoplanetary disks by analysing 
the SEDs of the young stars, taking advantage of the fact 
that the EXPORT optical and near-IR photometry were 
obtained simultaneously. This observational approach is 
appropriate since T Tauri and HAeBe stars vary markedly 
in these spectral regimes and a significant part of the to- 
tal luminosity of the object is radiated by the PMS stellar 
photosphere at these wavelengths. Among the stars in the 
EXPORT sample with measured Hipparcos parallaxes, the 
HAeBe stars HD 34282 and HD 141569 are the only ones 
whose positions fall below the zero-age main sequence in 
a logL — logTcff HR diagram or equivalent (e.g. van den 
Ancker et al. 1998, Weinberger et al. 2000). This result is 
difficult to reconcile with some observational results and 
with theoretical PMS evolutionary models (e.g. Yi et al. 
2001). 

In this paper we present an analysis of the stellar prop- 
erties of HD 34282 and HD 141569 and their circumstellar 
disks, based on EXPORT data complemented with new 
spectroscopic observations and data from the literature. 
The structure of the paper is as follows. In Section El we 
briefly review what is currently known about the stars. 
Sections El and 01 present the observations used in this 
work and their results. In Section [SI we do a comprehen- 
sive study of the stellar parameters and distances to the 
stars. In Sectional the disk models used to reproduce the 
observed SEDs are presented. In Section d we summarize 
the results of the work. 

2. The stars 

HD 34282: It is a HAeBe star with spectral type estimates 
in the range AO V - A3 V and a rotational velocity of 
V sin i = 129 ±8 km/s (Mora et al. 2001). It is variable with 
pubhshed my values between 9.8 and 10.11 mag (Sylvester 
et al. 1996; de Winter et al. 2001). Malfait et al. (1998) re- 
port a visual variability amplitude of ~ 2.5 mag, while an 



Algol-type minimum of Am ~0.8 mag is observed in the 
Hipparcos light curve, indicating that HD 34282 is most 
likely an UXOR-type object. Its SED shows a strong IR 
excess which is already noticeable at near-IR wavelengths 
(Sylvester et al. 1996; Malfait et al. 1998). The Hipparcos 
parallax is tt = 6.10 ± 1.63 mas, equivalent to a distance 
of 164^30 P*^- Based on this distance, van den Ancker et 
al. (1998) estimate a luminosity of 4.8 L©, a value which 
locates HD 34282 in an anomalous position in the HR di- 
agram, well below the expected luminosity of a star of its 
spectral type. Note that the error in the parallax, namely 
cr = 1.63 mas, is slighly larger than usual for Hipparcos 
objects, most probably due to the faintness of the star, af- 
fecting the estimate of the distance. A large keplerian disk 
around HD 34282 has been inferred from interferometric 
12 CO J = 2^1 line observations (Pietu at al. 2003). Those 
authors estimate a distance cZ = 400^}qo pc based on dy- 
namical considerations, which yields L„/Lq — 29^^^^ and 
M*/M0 = 2.1^0 2 for the star, much more consistent with 
the expected stellar properties of HD 34282. 

HD 141569: This star is considered to be a transi- 
tional object between the PMS HAeBe stars and the more 
evolved Vega-type systems. Optical and near-IR images 
reveal a circumstellar disk with a complex annular struc- 
ture extending up to a radial distance of about 600 AU 
(Weinberger et al. 1999; Augereau et al. 1999; Mouillet et 
al. 2001; Clampin et al. 2003). Mid-IR images trace the 
inner portion of the disk, up to about 100 AU (Fisher 
et al. 2000; Marsh et al. 2002). The disk in HD 141569 
shares many of the properties of dusty debris disks associ- 
ated with young MS stars, such as (3 Pic or HR 4796. On 
the other hand, the gas content of the disk is significantly 
higher than that of MS debris disks, as evidenced by the 
CO measurements at millimetre and near-IR wavelengths 
(Zuckerman et al. 1995; Brittain & Rettig 2002; Boccaletti 
et al. 2003). 

Concerning the star itself, HD 141569 has spectral type 
estimates between B9 and AOV and a very high rotational 
velocity, u sin i = 258±17 km/s (Andrillat et al. 1990; Mora 
et al. 2001 and references therein). The star has no photo- 
metric variability (Alvarez & Schuster 1981). Its spectrum 
exhibits double-peaked Ha emission superimposed on a 
broad and strong absorption line (Andrillat et al. 1990; 
Dunkin et al. 1997) as well as double-peaked emission of 
O I at 7772 and 8446 A (Andrillat et al. 1990); an outflow 
of material is suggested by the Na i D lines (Dunkin et 
al. 1997). Van den Ancker et al. (1998) estimate an age 
> 10^ years, while Weinberger et al. (2000) give an age 
of 5 Myr, on the basis of the stellar properties and age 
estimates of two nearby stars, which are most likely its 
companions and have T Tauri characteristics. Considering 
the Hipparcos distance, 99 pc, and the observed optical 
brightness and colours (e.g. de Winter et al. 2001), we de- 
duce an absolute magnitude of My = 1.8 mag, assuming 
an AOV spectral type and a typical interstellar extinction 
law (Rieke & Lebofsky 1985). This makes the star under- 
luminous, a property seemingly shared by the debris disk 
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Table 1. Telescopes, instruments, setups and observing modes. 



Telescope 


Instrument 


Range /Bands 


A/AA 


Observing mode 


WHT (4.2m) 


UES 


3800-5900 A 


~ 49000 


Echelle spectroscopy 


INT (2.5m) 


IDS 


5800-6800 A 


~5000 


Mid resolution spectroscopy 


NOT (2.5m) 


Turpol (photopol.) 


UBVRI 




Photopolarimetry 


GST (1.5m) 


Photometer+near-IR camera 


JHK 




Photometry 


CAHA 2.2m 


CAFOS (spectros.) 


3700-6440 A 


~2400 


Mid resolution spectroscopy 



Observing dates for the spectroscopic observations^ 





HD 34282 


HD 141569 


WHT/UES 


28, 29 January 99 


16, 17 May 98 (2+2) 

28, 29, 30 July 98 

28, 29, 30, 31 January 99 


INT/IDS 


24, 25, 26, 27, 28 October 98 


14, 15, 16, 17 May 98 




29, 30, 31 January 99 


29, 30, 31 July 98 
29, 30, 31 January 99 


CAHA 2.2/CAFOS 


3 February 03 


6 March 03 



f The observing logs for the photometric observations are given by Eiroa et al. 
(2001) and Oudmaijer et al. (2001). 



systems (3 Pic, HR 4796 and some other A-type young 
MS stars (Jura et al. 1998, Lowrance et al. 2000). We 
note, however, that HD 141569 has the largest IR excess 
among all these objects, it is the only object with emission 
lines in its optical spectrum, there is evidence of outflow- 
ing gas, its disk has a remarkable gas content, and it likely 
has the smaller evolutionary age. In addition, from a PMS 
evolution theory point of view, a 5 Myr old AO star {Tcs 
= 10000 K, see below) with solar metallicity would have 
an absolute magnitude My = 1.0-1.1 mag (e.g. Yi et al. 
2001), which is significantly brighter than the magnitude 
deduced from the Hipparcos parallax. 

3. Observations 

High resolution echelle (wavelength range 3800-5900 A) 
and intermediate resolution (wavelength range 5800-6700 
A) spectroscopy, optical photo-polarimetry and near- 
infrared photometry of HD 34282 and HD 141569 were 
obtained during the four EXPORT observational cam- 
paigns in May, July and October 1998 and January 1999 
carried out at the Canary Islands observatories (Eiroa 
et al. 2000). Details of the observations, instrumental se- 
tups and reduction procedures for each particular observ- 
ing mode are given in Eiroa et al. (2001), Mora et al. 
(2001) and Oudmaijer et al. (2001), so we do not repeat 
them here, though Table ^ provides an overview of the 
telescopes, instruments, setups and observing modes used 
during the EXPORT campaigns, plus the additional ob- 
servations carried out in 2003 (see below). The dates for 
the spectroscopic observations are also specified. 

The spectra of HD 34282 and HD 141569 are presented 
and analysed here for the first time, while the photometric 
data, though already published as such in the previous 
references, are collected here together and also discussed 
for the first time. We note that the data arc simultaneous 
or almost simultaneous, allowing a coherent study to be 
made of the photometric and spectroscopic variations, and 



the construction of the corresponding SEDs in the interval 
UBVRIJHK with contemporaneous data. 

In addition to the EXPORT observations, intermedi- 
ate resolution spectra in the range 3700-6200 A, with very 
high signal-to-noise ratio, have been taken for both stars 
with CAFOS (Calar Alto Faint Object Spectrograph) on 
the 2.2m telescope at Calar Alto observatory (Almerfa, 
Spain). The motivation was to obtain high quality Balmer 
line profiles, in particular H/3, H7 and H(5, in order to es- 
timate the stellar gravities (see Section O . The echelle 
spectra were inadequate for analysis of the Balmer line 
profiles: these are so wide that their wings extend up to 
three echelle spectral orders - the correction of the re- 
sponse function of the spectrograph and the splicing of 
adjacent orders are not accurate enough to consider the 
reconstructed profiles as reliable. CAFOS was equipped 
with a CCD SITe detector of 2048x2048 pixels (pixel size 
24 nm) and the grism Blue-100, centered at 4238 A, giv- 
ing a reciprocal linear dispersion of 88 A/mm (2 A/pixel). 
Three spectra, with exposure times 120, 1200 and 1200 s, 
were taken for HD 34282 on 5 February 2003, and two with 
exposure times 10 and 100 s were obtained for HD 141569 
on 6 March 2003. The usual bias, dark and dome flat- 
field frames were taken in each case. Standard procedures 
were used to process the data. The nominal wavelengths 
of the Balmer lines were used to self-calibrate the spectra 
in wavelength. 

4. Results 

4.1. Spectroscopy 

The high and intermediate resolution spectra of both stars 
show non- variable early-type photospheric spectra, inter- 
stellar (IS) narrow absorption lines and diffuse IS bands. 
In both cases, however, the most prominent features come 
from the circumstellar (CS) environment. Fig.^shows the 
average INT spectra of HD 34282 (January 99) and HD 
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Wavelength (A) 



Fig. 1. The normahzed average IDS/INT spectra of HD 
34282 and HD 141569. Some prominent lines are identi- 
fied. DIB stands for "diffuse interstellar band" . The spec- 
trum of HD 34282 is the average of the January 1999 run 
while the spectrum of HD 141569 is the average of the 
May 1998 run. 
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141569 (May 98) where some lines are identified. In the 
following, we describe the CS spectrum of both objects. 

HD 34282: The UES echelle spectra show a CS contri- 
bution in the Ca ii K line and in the Balmer lines H/3 and 
H7. The line He i at 5876 A is seen in absorption; it is 
unlikely to be of photospheric origin since HD 34282 has 
an ~A3 spectral type. Its strength and profile might show 
small variations (after a comparison of all the INT spec- 
tra), but this should be confirmed with higher resolution 
spectra. He i 5876 A has been observed with large vari- 
ations in HAeBe and T Tauri stars (e.g. Bohm & Catala 
1995, Johns & Basri 1995) and it has been explained in 
terms of magnetospheric accretion in the low-mass PMS 
stars (e.g. Oliveira et al. 2000), but in the case of HAeBe 
stars its origin is controversial (e.g. Bouret & Catala 2000). 

Ha always presents a variable double-peaked emission 
with a central absorption. Fig. El shows the observed Ha 
profiles after the subtraction of a spectrum of a standard of 
the same spectral type broadened to the rotational veloc- 
ity of HD 34282. The depth of the central absorption varies 
and its radial velocity, measured with respect to the stel- 
lar photosphere, shifts from an average value of 27.5±5.5 
km s~^ in October 98 to -63.3±14.6 km s~^ in January 
99. The radial velocities of the violet and red emission 
peaks and their separation also change from October 98 
to January 99 (the maximum peak separation is 319 km 
s^^ on 28 Oct 98 and the minimum is 117 km on 30 Jan 
99). The relative intensity of the emission peaks changes 
from day to day and a flip is observed from October 98, 
when the violet peak is stronger, to January 99, when the 
read peak is the stronger one. Furthermore, the blue wing 
shows a shoulder on 28 October 98 which could be due to 
a third velocity component. Finally, the total equivalent 



Fig. 2. Ha profiles for HD 34282 and HD 141569. Spectra 
of standard stars with spectral types A3 V and AO V, 
respectively, have been subtracted from the original spec- 
tra, therefore only the pure emission is shown in these 
panels. The lines have been shifted vertically by an ar- 
bitrary amount for clarity, so their wings should be put 
at the same level in order to compare the profiles. Each 
profile has been shifted horizontally 30 A with respect to 
the previous one as indicated in the abscissae. The vertical 
dotted lines mark the rest wavelength for Ha. The spectra 
have been ordered chronologically from left to right (see 
the dates of the observations in Table 

width varies from -3.0 to -5.7 A. Using the R band fiux 
from simultaneous optical photometry on 27 October 98 
and 30 January 99 we estimate a Ha line emission flux of 
7.4x10"^^ and 1.6x10"^^ erg cm~^ s^^, respectively. The 
variability in the emission flux and shape of Ha proves the 
changing conditions in the CS region where this line orig- 
inates. Grinin & Rostopchina (1996) studied a sample of 
HAeBe stars and found that double-peaked Ha emission 
likely arises in irregular gaseous circumstcUar disks rotat- 
ing close to the stars and seen nearly edge on. We notice 
that Pietu et al. (2003) find an inclination of the CO ke- 
plerian disk of i = 56° with respect to the plane of the 
sky. The change of the strength ratio between the violet 
and red peaks could reflect the relative contribution of the 
approaching and receding parts of the disk in the different 
spectra. More sophisticated accretion disk models have re- 
cently been developed by Tambovtseva & Grinin (2000) 
to explain the hydrogen lines of HAeBe stars in a context 
similar to the widely accepted magnetospheric accretion 
models for the low-mass T Tauri stars. 
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Object 


u 


B 


V 


R 


I 


J 


H 


K 


Refs. 


HD 34282 


10.47 


10.23 


10.03 


9.91 


9.76 


9.11 


8.40 


7.62 


1, October 98 




10.19 


10.03 


9.87 


9.77 


9.63 


9.15 


8.53 


7.84 


1, January 99 








10.11 






9.17 


8.32 


7.42 


2 




10.16 


10.05 


9.88 


9.79 


9.68 


9.02 


8.20 


7.44 


3 




10.15 


10.05 


9.89 


9.81 


9.71 








4 




10.11 


9.97 


9.81 


9.73 


9.62 








4 




10.51 


10.07 


9.94 


9.95 


9.79 








4 














8.99 


8.54 


8.08 


6 


HD 141569 


7.23 


7.20 


7.11 


7.01 


7.03 


6.70 


6.59 


6.53 


1 








7.09 






6.88 


6.86 


6.82 


2 








7.13 






6.87 


6.84 


6.80 


3 




7.24 


7.21 


7.10 


7.03 


6.96 








4 




7.28 


7.26 


7.18 
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(1) EXPORT, (2) Malfait et al. (1998), (3) Sylvester et al. (1996), (4) de Winter et al. (2001), 
(5) Penprase (1992), (6) Coulson et al. (1998) 
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Fig. 3. Details of the Ha profiles for HD 34282 and HD 
141569. The most extreme cases have been taken for each 
star to illustrate the complex behaviour of this line. The 
variations in the absorption profiles appear in both wings 
but are more pronounced in the blue one. 

HD 141569: Double-peaked Ha emission is observed in 
the spectrum of this star. Such double-peaked emission, 
though much less pronounced, is also observed superim- 
posed to the H/3 photospheric line in all UES spectra. The 
heliocentric radial velocity of the star obtained from the 
UES spectra is -13.1 km s~^, which is in good agreement 
with previous estimates (Dunkin ct al. 1997, Brittain et 
al. 2003), given the velocity resolution of the UES spec- 
tra, 6 km s^^. The violet and red Ha emission peaks have 
— 113.2±3.9 km/s and 128.3±5.0 km/s, respectively, while 
the central absorption lies at a velocity of 11.9±4.2 (aver- 



age values and rms error of the 10 INT spectra, velocities 
measured with respect to the stellar photosphere). Dunkin 
et al. (1997) find agreement between the stellar and the 
Ha central absorption velocities; the discrepancy in our 
values must be treated with caution given the spectral 
resolution of the individual INT spectra, ~ 50 km s~^. 
The observed velocity separation between both emission 
peaks is in good agreement with the results of Dunkin et 
al. (1997). The emission profiles are quite similar in all 
spectra, although they vary noticeably. The violet peak 
is always stronger than the red one, but their relative in- 
tensity changes (Fig. |2Il; there is also a variation in the 
Ha blue wing (Fig. |3Jl. In addition, the equivalent width 
varies, ranging from -6.14±0.03 A in July 1998, when the 
violet peak was larger, to -5.63±0.25 A in January 1999, 
when the violet peak was smaller. The average Ha line 
flux is 2.1 X 10^^^ erg cm~^ s^^. 

All INT spectra show [O i] in emission, relatively 
strong at 6300 A and very weak at 6364 A. The emission 
is constant and has an equivalent width W([0 i] 6300 A) 
0.13 A; its flux is 4.6 x 10~^^ erg cm"^ s"^ The peak of 
the emission occurs at the same radial velocity as the Ha 
absorption peak. [O i] emission is usually interpreted as 
an indicator of low density winds in HAeBe stars (Bohm 
k Catala 1994; Corcoran & Ray 1997). 

4.2. Photometry and the spectral energy distribution 

EXPORT optical UBVRI and near-IR JHK magnitudes 
of HD 32482 and HD 141569 are given in Oudmaijer et al. 
(2001) and Eiroa et al. (2001). HD 34282 varies in both 
wavelength regimes from October 98 to January 99, but no 
variability is detected on time scales of days. HD 141569 
remains constant on all time scales within the photometric 
uncertainties (typically a few percents in any photometry 
band). Table [5] shows the mean magnitudes of HD 34282 
as observed in October 98 and January 99, as well as the 
mean magnitudes of HD 141569. For the sake of com- 
parison, photometric results taken from the literature are 
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Table 3. Absolute dereddened fluxes for the SEDs 







HD 34282 


HD 141569 




Band 


log A 


log Fa 


log Fa 






(/im) 


(W m-2/xm-i) 


(W m^^^m" 


) 


WE 


-0.86 


- 


-9.90±0.02 


(4) 


lUE 


-0.75 


— 


-9.79±0.02 


(4) 


WE 


-0.66 


— 


-9.94±0.02 


(4) 


WE 


-0.59 


- 


-10.01±0.02 


(4) 


WE 


-0.53 


- 


-10.08±0.02 


(4) 


Wal W 


-0.49 


-11.45±0.02 (5) 


— 




Wal U 


-0.44 


-11.44±0.02 (5) 


— 




Wal L 


-0.42 


-11.23±0.02 (5) 


— 




Wal B 


-0.37 


-11.11±0.02 (5) 


— 




Wal V 


-0.26 


-11.32±0.02 (5) 


- 




U 


-0.44 


-11.35±0.01 (1) 


-10.02±0.02 


(1) 


B 


-0.36 


-11.12±0.01 (1) 


-9.88±0.02 


(1) 


V 


-0.26 


-11.31±0.01 (1) 


-10.12±0.02 


(1) 


R 


-0.15 


-11.62±0.01 (1) 


-10.45±0.02 


(1) 


I 


-0.01 


-11.90±0.01 (1) 


-10.82±0.04 


(1) 


J 


0.10 


-12.15±0.02 (1) 


-11.22±0.02 


(2) 


H 


0.21 


-12.32±0.02 (1) 


-11.64±0.02 


(2) 


K 


0.35 


-12.54±0.02 (1) 


-12.12±0.02 


(2) 


ISO 


0.39 


- 


-12.27±0.02 


(6) 


ISO 


0.51 


- 


-12.66±0.02 


(6) 


L 


0.58 


-12.75±0.02 (3) 


-12.92±0.02 


(3) 


ISO 


0.63 




-13.14±0.02 


(6) 


M 


0.68 


-13.12±0.02 (3) 


-13.25±0.02 


(3) 


ISO 


0.77 


— 


-13.46±0.02 


(6) 


ISO 


0.83 


— 


-13.54±0.02 


(6) 


ISO 


0.89 


— 


-13.39±0.02 


(6) 


ISO 


0.94 


— 


-13.61±0.02 


(6) 


ISO 


0.99 


— 


-13.78±0.02 


(6) 


ISO 


1.03 


- 


-13.82±0.02 


(6) 


ISO 


1.06 


- 


-13.84±0.02 


(6) 


IRAS 


1.08 


-13.82±0.10 (7) 


-13.92±0.05 


(7) 


IRAS 


1.40 


-14.09±0.10 (7) 


-14.03±0.07 


(7) 


IRAS 


1.78 


-14.03±0.10 (7) 


-14.32±0.09 


(7) 


IRAS 


2.00 


-14.50±0.10 (7) 


-15.00±0.10 


(7) 


0.45 mm 


2.65 


-16.72±0.23 (3) 


- 




0.8 mm 


2.90 


-17.72±0.07 (3) 






1.1 mm 


3.04 


-18.34±0.09 (3) 


< -19.05 (3) 




1.3 mm 


3.12 


-18.71±0.09 (8) 






1.35 mm 


3.13 




-20.05±0.02 


(9) 


2.6 mm 


3.41 


-19.99±0.13 (10) 






3.4 mm 


3.53 


-20.89±0.06 (8) 







(1) EXPORT, (2) Malfait et al. (1998), (3) Sylvester et al. 
(1996), (4) INES Database, (5) de Geus et al. (1990), (6) 
ISO PHOT-S, TDT No. 62701662, (7) IRASVSC (colour cor- 
rected), (8) Pietu et al. (2003), (9) Sylvester et al. (2001), (10) 
Mannings & Sargent (2000) 



also shown. In general, typical intrinsic error bars in the 
EXPORT UBVRI magnitudes are of the order of ±0.04 
mag and always less than 0.10 mag; for JHK they are less 
than 0.05 mag (see Eiroa et al. 2001 and Oudmaijer et al. 
2001 for specific values). 

HD 34282 was brighter in the optical in January 99 
than in October 98 (the values of V are 10.02, 10.03 
and 9.87 for 24, 28 October 98 and 30 January 99, re- 
spectively) . The January 99 magnitudes likely correspond 



to the maximum optical brightness level of this star, as 
comparison with other results (Table |2) and with the 
Hipparcos light curve suggests. The near-IR brightness be- 
haves in the opposite way: the star was brighter in October 
98. The optical variability amplitude decreases with wave- 
length, as expected from the usual mechanisms invoked 
to explain the photometric variability in PMS stars (e.g. 
Herbst & Shevchenko 1999), but in the near-IR the max- 
imum variability amplitude occurs in K. The cause of 
this different behaviour might be due to the fact that the 
optical variability is due to changes affecting the stellar 
photosphere while the near-IR fluxes reflect the contribu- 
tion of the star and the circumstellar disk. A simultane- 
ous ant icor relation of the variability in both wavelength 
regimes has also been found in other pre-main sequence 
stars, though on a much shorter time scale (Eiroa et al. 
2002; Grinin et al. 2002). 

Concerning HD 141569, our optical photometry and 
previous results are in excellent agreement. The apparent 
discrepancy between our near-IR magnitudes and those 
obtained by Sylvester et al. (1996) and Malfait et al. 
(1998) is likely due to the fact that the companion star 
HD 141569B was at least partly in our beam (Weinberger 
et al. 2000). Once our measurements are corrected for this, 
the agreement is satisfactory. 

For both stars EXPORT photometry together with 
data from near ultraviolet to millimetre wavelengths have 
been used to delineate the spectral energy distributions, 
which are, to our knowledge, the most complete ones 
to date. Figs. El and EH show the observed spectral en- 
ergy distributions of both stars. Table |21 gives dereddened 
fluxes for each star and corresponding references (they 
are also plotted in Figs. 151 and llOfl. E{B — V) values of 
0.05 for HD 34282 and 0.12 for HD 141569 (see sub- 
section 15.11 for the origin of these values) were used for 
the reddening corrections. The extintion curve of Rieke 
& Lebofsky (1985), from 4400 A to longer wavelengths, 
and that of Steenman & The (1991) for shorter wave- 
lengths were used. We have used the UBVRI (one set 
of values) and the mean JHK (three sets) photometry 
from January 99 for HD 34282 because they likely repre- 
sent better the stellar photosphere, i.e. when the star was 
brighter in the optical. For HD 141569, we used the mean 
EXPORT [75 ra/ photometry and JiTTf magnitudes from 
Malfait et al. (1998). UBVRIJHK magmt-ades have been 
converted into fluxes using the absolute calibration given 
by Schmidt-Kaler (1982). lUE spectra in the range 1250- 
3000 A of HD 141569 are taken from the INES database^ 
(Solano et al. 2001). No lUE observations are available 
for HD 34282. Walraven photometry for HD 34282 was 
not contemporaneous with EXPORT observations and has 
only been used in subsection l6.2.ll for a determination of 
the disk mass accretion rate. In the infrared, we have used 
colour-corrected IRAS fluxes and ISO scientifically vali- 
dated 'Off-Line Processing' products (OLP) taken from 
the Data Archive at the ISO Data Centre. Submillimetre 

^ ,http : // ines . laef f . esa.es/ ^ 
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and millimetre continuum fluxes have been taken from the 
references listed in Tabled the inclusion of the continuum 
fluxes at these long wavelengths alleviates, to some extent, 
the uncertainties in the flux calibration of the infrared ISO 
pipeline-reduced fluxes. 

5. Stellar parameters of the stars and their 
distances 

As we mentioned above, both stars have been reported to 
have anomalous positions in the HR diagram, falling be- 
low the main sequence. However, the observational data 
show that HD 34282 is a highly variable HAeBe star with 
an active CS environment. Concerning HD 141569, though 
its disk shares similarities with some debris disk systems, 
which are found to be located below the main sequence, 
the star does show observational properties which suggest 
an earlier evolutionary stage. In particular, our data show 
that the object retains a non-negligible CS activity, as 
evidenced by the [O i] emission and the variable Ha emis- 
sion. Taking into account these considerations, both stars 
should be located above or on the main sequence, in bet- 
ter agreement with normal expectations from their photo- 
metric and spectroscopic properties, and stellar evolution 
theory. In this section we carry out a detailed study of the 
physical parameters of the stars leading us to re-consider 
their location in the HR diagram. The determination of 
the effective temperature, gravity and stellar metallicity 
will allow us, by using the appropriate set of evolutionary 
tracks and isochrones, to estimate the luminosity, mass 
and age of the stars, as well as their distances by means 
of the measured stellar fluxes. 

5.1. Stellar parameters: effective temperature, gravity, 
metallicity and rotational velocity 

Kurucz (1993) model atmospheres are used to estimate 
the physical parameters describing the photospheres of 
HD 34282 and HD 141569. We use the following iterative 
procedure. First, an initial value of the stellar gravity is 
obtained from the comparison of synthetic profiles, com- 
puted with solar metallicity, with the observed H/?, H7 and 
H5 profiles extracted from the CAFOS observations. For 
this purpose, a tentative value of Toff, based on the spec- 
tral type of each star, is taken. The gravity is obtained 
by measuring the width of the observed Balmer profiles 
at one intensity level below the normalized continuum, 
namely I = 0.80; this has been chosen to avoid potential 
artifacts in the wings near the continuum arising from the 
normalization of the original spectra. The widths of the 
Balmer lines are then compared with those measured on a 
grid of synthetic profiles computed for different gravities^. 

^ The convolution of the theoretical profiles with the instru- 
mental response changes the results by a minute quantity. For 
HD 141569, which is the star with less intrinsic uncertainties 
in the determination of the gravity, the relative difference with 
the value in TableE]is about 0.5%. 
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Fig. 4. The Balmer lines H/3, H7 and H(5 for HD 34282 and 
HD 141569 (solid lines) extracted from the spectra taken 
with CAFOS on the 2.2m telescope at CAHA. The dotted 
lines are the Kurucz synthetic profiles corresponding to 
the Tcff, \ogg* and metallicity listed in Table 01 for each 
star. The discrepancies in the higher parts of the wings 
arise from contamination with weak metalic lines and also 
from the difficulties in normalizing the original spectra. 

A linear interpolation is done taking the widths of each 
Balmer line as independent variables and the gravity as 
the unknown. The result for log 17, is the mean value of 
all the interpolations with an uncertainty corresponding 
to the standard deviation. 

With this value of logg» fixed, we re-estimate the efec- 
tive temperatures of the stars and their mctallicities. This 
is done by means of a systematic comparison of synthetic 
line profiles, obtained using Kurucz 's models of different 
temperatures and mctallicities, with the high resolution 
WHT/UES spectra. The regions of effective temperatures 
from 8500 to 9750 K for HD 34282, and from 9750 to 10500 
K for HD 141569 were explored, which approximately en- 
compass the values corresponding to the spectral types 
reported for both stars. For each temperature, different 
mctallicities were used to synthesize spectral lines in sev- 
eral selected spectral intervals. The synthetic spectra were 
rotationally broadened with the corresponding v sin i val- 
ues for each star (Mora et al. 2001). A microturbulence 
velocity of 2 km/s was used throughout. 

For HD 34282 the lines chosen are: 4404.75, 4957.60 
A (Fe i), 4508.29, 4515.34, 4583.84, 4629.34 A (Fe 11) and 
4501.27 A (Ti 11). For HD 141569 the projected rotational 
velocity is so high that it is difficult to find individual 
metalhc lines; instead the regions 4260-4280 A, 4375-4400 
A and 4540-4560 A, which contain a large number of lines, 
have been selected. For both stars the lines 3933.66 A (Ca 
11), 4481.13 A(Mg 11), and the Si 11 doublet at 6347.11, 
6371.37 A have also been used (the Si 11 doublet is clearly 
observed in the INT spectra). It is found, particularly in 
the case of HD 34282, that different (Teg, [Fc/H]) pairs 
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fit the observed spectra reasonably welL The comparison 
shows that higher effective temperatures require higher 
metallicities to fit a given spectral line. This (Toff, [Fe/H]) 
degeneracy can be solved by using the observed SEDs, in 
particular the t/iJViZ/ photometry. 

For each (Tea, [Fe/H]) pair, which provides an ap- 
parent satisfactory fit to the observed spectra, we have 
computed the corresponding synthetic Kurucz SED in the 
spectral interval from 3000 to 10000 A, i.e., the range from 
U to /. The observed fZB V7?/ fluxes from each star have 
been dereddened at small steps of E{B—V) and compared 
to the Kurucz synthetic models (we have covered the in- 
terval 0.0 < E{B-V) < 0.30 with a stepsize of 0.01). 
The observed and synthetic photometry were normalized 
to R. It turns out that the overall SED shape is signifi- 
cantly more sensitive to the effective temperature than to 
the metallicity; thus, this method is particularly useful to 
pinpoint the best Tee, and consequently the (Toff, [Fe/H]) 
solution. The best Teff value fitting the data is chosen by 
means of a minimum criteria. Once the best values for 
Teff and [Fe/H] have been found, the value of the gravity 
is re-calculated using the procedure described above with 
synthetic Balmer profiles corresponding to this metallicity. 
The process is iterated until convergence is achieved. 

The optimum solutions found are (see also Table ^ : 
HD 34282, Teff = 8625 K, [Fe/H] =-0.8, logy* = 4.20; HD 
141569, Teff = 10000 K, [Fe/H] =-0.5, logg* = 4.28. The 
corresponding spectral types for these temperatures and 
gravities are A3 V and B9.5 V respectively. The relevant 
result is that both stars have low metal abundances. In 
addition, projected rotational velocities, wsini, of 110 km 
s-i for HD 34282 and 236 km s^^ for HD 141569 provide 
a better agreement with the observed profiles than the 
vs'mi values reported by Mora et al. (2001). Finally, the 
best fit gives E{B -V) = 0.05 mag for HD 34282 and 
0.12 mag for HD 141569. We note that Mora et al. (2001) 
deduced E{B — V) = 0.15 ± 0.10 mag for the latter star 
from the graphite interstellar absorption band at 2200 A 
observed with the lUE, also Oudmaijer et al. (2001) found 
Ay = 0.34 ± 0.20, which implies E{B-V) = 0.11 ± 0.06; 
both values of the extinction are in very good agreement 
with our present estimate. Figs. |5l and |H1 show the compar- 
ison between the observed and best synthetic spectra for 
several spectral lines and ranges; synthetic spectra using 
solar abundances are also plotted for comparison. Both fig- 
ures clearly show that the observed spectra of HD 34282 
and HD 141569 are much better reproduced with our low 
metallicity estimates and, therefore, a low metal content 
in the photospheres of both stars is strongly suggested. We 
note that HD 34282 has some similarities with HD 37411, 
another low metaUicity PMS star (Gray & Corbally 1998). 
With respect to HD 141569, Gray & Corbally (1998) in- 
dicate that it is a 'mildy metal-weak' star, in agreement 
with our results, though Dunkin et al. (1997) give a solar 
abundance. 



Table 4. Stellar parameters determined. 

HD 34282 

Effective temperature 8625±200 K 

Spectral type A3 V 

wsini 110±10km/s 
Metallicity [Fe/H] = - . 8 ± . 1 

Log fl. 4.20±0.20 



E{B - V) 




0.05 


Stellar flux (K) 
IR Excess (Fir) 
Air, 


(3.61 ± 0.44) X 10"'^ W m"^ 
(2.28 ± 0.40) X 10~^^ W 
1.25 /im 




[±Alog<7.] 


[iATefi'] 


L/Lq 

Mass (M0) 
Age (Myr) 
Distance (pc) 


13.64]:--2 
1 50-0.07 

'-'•^-L-2.58 
348+129 


13.64+i-^« 

1 50+0.04 
J^->^^~0.05 
541-0.41 
"•^J^+0.46 

348t^^ 


HD 141569 



Effective temperature 10000±200 K 

Spectral type B9.5 V 

vsini 236±15km/s 

MetaUicity [Fe/H] = - . 5 ± . 1 

Log g» 4.28±0.04 



E{B - V) 




0.12 


Stellar flux (F,) 


(7.16 ±0.24) X 10-"" W m"^ 


IR Excess (Fir) 


(1.27 ±0.27) X 10'"2 W m-2 


Air 




3.25 /im 




[±Alog<7, 


] [iATefi-] 


L/Lq 

Mass (M0) 
Age (Myr) 
Distance (pc) 


25.77;^:^^ 
2.00«;gl 

4.7i«:?^ 

108;^ 


' ' -2.54 

2 00+°"'' 
^■'J'-'_o.o5 

A 71-0.31 

*-'-'^+0.33 

1081^ 



5.2. Stellar fluxes and infrared excesses. Luminosities, 
masses, ages and distances. 

Once the stellar parameters defining the photospheres of 
HD 34282 and HD 141569 have been determined, we can 
calculate the stellar fluxes and infrared excesses. In addi- 
tion, we can locate the stars in a HR diagram using appro- 
priate evolutionary tracks and estimate their distances. 

The stellar flux, F^,, is computed by integrating the 
best-fit Kurucz model normalized to the dereddened ob- 
served i?-band flux, as given in Table |31 The infrared ex- 
cess, Tift, is calculated by subtracting the flux under the 
tail of the best-fit Kurucz model from the observed fluxes, 
and integrating the result from Xib,, the wavelength where 
the SED separates from the tail of the Kurucz model (see 
Figs. and I10|l . Extrapolations to long wavelengths fol- 
lowing the formula by Chavarria et al. (1981) have been 
included in the integration. Table 0] shows the integrated 
fluxes for the stars and their IR excesses. HD 34282 has a 
Pm/F^ ratio of 0.63 ± 0.19, which suggests that some in- 
trinsic luminosity might be arising from the disk; Tir/T» 
values of ^^0.25 and ~0.5 are expected from flat and flared 
optically thick reprocessing disks, respectively (Kenyon & 
Hartmann 1987). In the case of HD 141569, we find Tjr/T, 
= 0.018 ± 0.004 and a passive disk with relatively low op- 
tical depth is suggested. This figure is one order of magni- 
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Fig. 5. Comparison of observed and synthetic spectral features of HD 34282. From top to bottom and left to right the 
following lines have been plotted: Ca ii K (3933.66 A), Mg ii (4481.13 A), Si ii (6347.11, 6371.37 A), Fe ii (4583.84 
A), Fe I (4957.60 A) and Fe ii (4629.34 A). Observed spectra are plotted as solid lines. Synthetic spectra computed 
with Toff = 8625 K, [Fe/H]=-0.8 and \ogg^= 4.20 and broadened with t;sini = 110 km/s are plotted as dashed lines. 
Dotted lines represent synthetic spectra computed with the same Tgff, logg, and vsini but using solar abundances. 
In the lower panels, the synthetic unbroadened spectra are shown to illustrate the lines giving rise to the observed 
profiles. 



tudc higher than typical values in Vega-type stars, (10~^ 
to lO^'', Sylvester et al. 1996). We note that the ratio 
Fik/F^ — 8 X 10~^ given by Zuckerman et al. (1995) and 
usually quoted in the literature does not include the 1.35 
mm flux measured by Sylvester et al. (2001). If we remove 
that flux from our estimates, we recover the figure given 
by Zuckerman et al. 

Knowledge of the stellar metallicities is relevant 
since changes in the metal abundance, Z, produce non- 
negligible changes in the relative positions of the evolu- 
tionary tracks. In particular, a decrease in the metallicity 
moves one given track with a fixed value of the stellar 
mass, to lower luminosities, higher gravities and higher 
temperatures, apart from changing its overall shape. This 
effect is illustrated in Fig. [T] for the log 5^ - log Toff HI^ 
diagram. The left two panels of Fig. [7| show the positions 
of the stars in the HR diagram according to the results 
found for Tcff and log^.^. The evolutionary tracks for a 
scaled solar mixture from the Yonsei-Yale group (Yi et 
al. 2001) - Y^ in their notation - have been used in this 
work. The Y^ tracks with Z = 0.004 and Z = 0.007 are ap- 
propriate for HD 34282 and HD 141569, respectively. For 



comparison, the tracks with solar metallicity {Z = 0.02) 
are also included in the plot'^. 

The corresponding values of the effective temperature 
and gravity are compared with the set of tracks and values 
for the luminosity, mass and age are derived. The position 
of HD 34282 and HD 141569 in the HR diagram log L/Lq 
- logTeff are plotted in the two panels at the right hand 
side of Fig. \7\ The luminosities, combined with the to- 
tal stellar fluxes obtained, imply distances of 348 pc 
for HD 34282, and 108 pc for HD 141569. All the results 
are summarized in Table ^ The results for the luminos- 
ity, mass, age and distance are given in two columns for 
each star: in the first one the error bars derive from the 



In the original set, there are no tracks for Z — 0.003 and 
0.006, the appropriate values for each star considering that 
the relationship Fe/Fe0=X/XQ holds (in number if particles), 
where X means any element with atomic number larger than 
3. The closest sets of tracks are those with Z = 0.004 and 0.007 
respectively. We have preferred to use original tracks kindly 
provided by the authors than interpolating between adjacent 
grids to build new sets of tracks and isochrones. In any case 
the differences are expected to be very small. 
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Fig. 6. Comparison of observed and synthetic spectral features of HD 141569. In the top panels, the spectral lines are 
the same as in Fig.|31 In the bottom panels, three spectral regions, instead of individual lines, are displayed. Observed 
spectra are plotted as solid hues. Synthetic spectra computed with reff = 10000 K, [Fe/H]=-0.5 and logg*= 4.28 and 
broadened with ?;sini = 236 km/s are plotted as dashed lines. Dotted lines represent synthetic spectra computed with 
the same Teff, \o'gg^ and wsin? but using solar abundances. In the lower panels, the synthetic unbroadened spectra 
are shown to illustrate the lines giving rise to the observed profiles. 



uncertainties in loggi*, keeping Toff fixed; in the second, 
the error bars derive from the uncertainties in Tes, keep- 
ing log fixed. The uncertainties for HD 34282 are larger 
than those for HD 141569 because of the larger uncertain- 
ties in the determination of the gravity. The vertices of 
the 'diamond-like' boxes appearing in the two right-hand 
panels of Fig. [7| correspond to the tips of the error bars in 
the panels on the left. 

Fig. [3 shows that both stars lie above the main se- 
quence in positions in agreement with their observational 
properties and current PMS evolution theory. The anoma- 
lous positions in the HR diagram previously reported for 
these stars referred to standard sets of tracks or photo- 
metric parameters computed with solar metallicity and a 
very doubtful distance in the case of HD 34282. For this 
star the anomalous position in the HR diagram is resolved 
both by the low metallicity and an appropriate distance 
estimate. In the case of HD 141569, its anomalous posi- 
tion originated only from the use of the wrong metallicity 
in the definition of the main sequence^. 

* Note that in the two left-hand-side diagrams of Fig.|7|both 
stars also occupy correct positions in the HR diagram when 
compared with tracks with solar metallicity. If these tracks 



The distance we obtain for HD 34282 (348l^f pc) is 
comparable to that given by Pietu et al. (2003) {•iOOtHo 
pc). Both values are remarkably similar, and are more 
than 2(7 larger than the Hipparcos value. This leads us 
to consider that the Hipparcos parallax for this object 
has been overestimated. The stellar parameters quoted by 
Pietu et al. (2003) are compatible with ours, although we 
note their larger errors and their use of indirect methods 
to estimate the temperature and spectral type of the star; 
in addition, they use solar metallicity tracks which are not 
appropriate for this star. In the case of HD 141569 our dis- 
tance is comparable to the Hipparcos one (99l^ pc). We 
note that the age estimate by Weinberger et al. (2000) for 
the two nearby low mass stars is similar to that deduced 
by us for HD 141569 directly. This result reinforces the 
idea that the three stars are real companions. 

Finally, estimates of the extinctions and distances can 
also be made using the absolute magnitudes and colours 
provided by the isochrones. Taking the values of mass, 
age and metallicity given in Table 01 and the observed val- 
ues of B and V used to construct the SEDs (see subsection 

were used to obtain values of masses, luminosities and ages, 
these quantities would be incorrect. 
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Fig. 7. Positions of the stars in the HR diagram according to the results shown in Table^l Left. Upper panel: Position 
of HD 34282 in the HR diagram log 5* - log Toff. The PMS tracks with Z = 0.004 (solid hnes) and Z = 0.G2 (solar, 
dotted lines) are plotted. Lower panel: same as in the upper panel for HD 141569. The tracks correspond to Z = 0.007 
(solid lines) and Z — 0.02 (dotted lines). Right. Upper panel: Position of HD 34282 in the HR diagram \ogL/LQ - 
logTcff. The PMS tracks are those for Z = 0.004 (solid lines). Lower panel: same as in the upper panel for HD 141569. 
The tracks correspond to Z ~ 0.007 (solid lines). In both cases we have labelled the tracks with the corresponding 
masses and the isochrones for 4, 6, 8 and 10 Myr are shown with dotted lines. 



14.21 and Table [21, the extinctions E{B — V) obtained are 
0.04 for HD 34282 and 0.06 for HD 141569. The former 
value especially is in excellent agreement with that given 
by the analysis. The absolute magnitudes My provided 
by the tracks are 2.04 and 1.37 respectively; these, com- 
bined with the observed values of V corrected with the 
above extinctions, yield distances of 347 pc for HD 34282 
and 129 pc for HD 141569, again in good agreement with 
the results shown in Tabled 

6. Modelling the SEDs: the irradiated accretion 
disk models 

The origin of the IR excess in HAeBe stars has been a 
matter of intense debate (see e.g. Waters & Waelkens 
1998). Some authors argued in favor of models with an ap- 
proximately spherical dusty envelope of low optical depth 
(Berrilh et al. 1992; di Francesco et al. 1994; Pezzuto 
et al. 1997; Miroshnichenko et al. 1997), others consid- 



ered circumstellar disks (Hillenbrand et al. 1992; Chiang 
et al. 2001; Natta et al. 2001) or proposed the existence 
of both (Natta et al. 1993; Miroshnichenko et al. 1999). 
Ultimately the existence of disks around HAeBe stars was 
firmly established by direct imaging at millimetre wave- 
lengths (Mannings et al. 1997; Mannings & Sargent 1997, 
2000; Testi et al. 2001; Pietu et al. 2003), in the optical 
(Grady et al. 1999, 2000) and in the IR (Close et al. 1997; 
Jayawardhana et al. 1998; Weinberger et al. 1999). Some 
authors have been able to fit the complete SEDs of HAeBe 
stars with passive irradiated disk models from Chiang & 
Goldreich (Chiang et al. 2001; Natta et al. 2001). 

6.1. The models 

In this work we use the self-consistent irradiated accretion 
disk models from D'Alessio et al. (1998, 1999 and 2001) 
to fit the SEDs of HD 34282 and HD 141569. The meth- 
ods used here were applied succesfully to fit the SEDs of 
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young T Tauri stars with an approximate age of 1 Myr. 
This approach is vahd for these two stars since they are 
young (see Section I^J and at least one of them seems to 
be still actively accreting. Since the calculations yield the 
vertical structure and emission properties of the disk self- 
consistently with the stellar parameters, and we have char- 
acterized the stellar photospheres with precision, a good 
fit to the SED would provide a physically based picture 
of the system without the use of ad-hoc parametrizations 
for the disk temperature or surface density profiles. 

The models involve the following assumptions: the disk 
is in a steady state (M= dAI/dt is constant), it is geomet- 
rically thin {H/ R<<1, where H is the scale height of the 
disk and R is the radial distance), the accretion viscosity 
is computed with v = a H Cg, following the a-prescription 
from Shakura & Sunyaev (1973), in which Cg is the sound 
speed. The dust and gas are well mixed in the whole disk 
with the usual dust to gas mass ratio of 1/100. For the 
dust we use a grain size distribution n{a) — uq a^^ with 
p values of 2.5 or 3.5 and a minimum grain size of 0.005 
/zm. Possible values for the maximum grain size are 1 /im 
to 10 cm. The abundances are those given by Pollack et al. 
(1994). The radiation field is considered in two separate 
regimes (one characteristic of the disk local temperature 
and one characteristic of the stellar effective temperature) 
as in Calvet et al. (1991, 1992) and finally the radiative 
transfer is done by solving the first two moments of the ra- 
diative transfer equation with the Eddington approxima- 
tion (see D'Alessio et al. 1998, 1999 and 2001 for further 
details). 

The disk is the solution of the detailed vertical struc- 
ture equations for each annulus, and the synthetic SEDs 
are obtained by solving the radiative transport equation 
in rays parallel to the line of sight and integrating (in solid 
angle) the disk's emission. 

Note that, although the models we have applied are 
amongst the most sophisticated available, they have lim- 
itations inherent in the complexity of the physics and 
the computational techniques. For example, they assume 
a unique dust size distribution with a maximum grain 
size flmax and the same viscosity parameter a throughout, 
which in some real cases precludes fitting the observed 
SED with a single component. For HD 34282 a combina- 
tion of two disk models is needed: one for the long wave- 
length part of the SED, with high viscosity and a maxi- 
mum grain size of 1 cm (subsection Ifi. 2.211 and a second 
one for the mid- infrared region of the SED, with lower 
viscosity and a maximum grain size of 1 /im (subsection 
I6.2.3() : the first one will be referred to as the mid-plane 
disk and the latter as the surface disk according to their 
vertical scale heights, of which the first has a lower value. 
We found that these two models were the only combina- 
tion able to reproduce the whole SED of this star, implying 
that some vertical grain segregation may be present in the 
disk around HD 34282. In addition, to explain the bump 
at 3 /im we propose the presence of a wall-like structure 
at the inner edge of the disk (subsection I6.2.4|l . For HD 



141569 a single disk model accounts for the whole SED 
('subsection l6.3|l . 

6.2. The disk model for HD 34282 

The disk around HD 34282 shows signs of activity: opti- 
cal and near-IR variability. Ha variable emission and high 
fractional IR excess. Hence, to fit the SED of HD 34282, 
we assume that the disk is heated by irradiation from the 
central star and viscous dissipation from the mass accre- 
tion onto the star. 

For the central star the following stellar parameters are 
assumed (TableEJ: Teff = 8625 K, i?* = 1.66 Rq, M* = 1.59 
Mq, = 13.64 L0 and d = 348 pc, where the stellar ra- 
dius was calculated from the stellar mass and the surface 
gravity. For the disk we take an inclination angle ofi — 56° 
(Pietu et al. 2003) and an outer radius of 705 AU, which 
is the disk radius given by these authors scaled from their 
distance of 400 pc to our value of 348 pc. The inner disk 
radius is calculated self-consistently using the position of 
the near-IR bump, the result being 0.31 AU (see subsec- 
tion EHiSIl • 

To find the remaining input disk parameters we need 
to know the mass accretion rate towards the star (which 
will be derived from ultraviolet Walraven photometry), 
the disk viscosity a (which will be fitted to raise the ap- 
propriate flux at millimetre wavelengths once we know the 
mass accretion rate) and the dust properties (related to 
the slope in the millimetre range and the IR excess shape 
at shorter wavelengths). 

6.2.1. The mass accretion rate in the disk of HD 34282 

The fractional IR excess, Fir/F*, of HD 34282, namely 
0.63 ± 0.19 (Table EI), is a little bit larger than the val- 
ues measured in the SEDs of the HAeBe stars studied by 
Natta et al. (2001), all of them being smaller than 0.45. 
These authors reproduce the SEDs assuming that no ac- 
cretion is present, the reprocessing of radiation being the 
mechanism invoked to explain the IR excess. We propose 
that, given the more prominent IR excess of this star com- 
pared to those of the HAeBe stars in Natta et al. (2001) 
and the higher fractional IR excess (slightly larger than 
0.5), accretion may be present in this object. For this rea- 
son, we first try to estimate the mass accretion rate to- 
wards HD 34282. 

This parameter has been studied using Walraven ul- 
traviolet photometry of HD 34282 (de Geus et al. 1990) 
and shock emission models by MuzeroUe et al. (2003a), 
who apply to HAeBe stars the magnetospheric accretion 
scenario from Calvet & GuUbring (1998) and GuUbring et 
al. (2000). This theory explains the UV excess emission 
seen in many T Tauri stars as arising from an accretion 
shock. In HAeBe stars, the emission from the shocked at- 
mosphere modifies the underlying photosphcric emission 
reducing the Balmer jump with respect to the normal 
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Fig. 8. Best-fit shock model fitting to the Walraven pho- 
tometry of HD 34282. The open and filled circles are 
observed and dereddened Walraven photometry, respec- 
tively. The open squares are synthetic photometry com- 
puted with the total (stellar plus shock) flux. The solid 
line is the Kurucz model plus the accretion shock flux for 
the accretion column energy flux F and filling factors / 
specified in the plot. Shifted a small amount below this 
line would lie the Kurucz model for the stellar parame- 
ters of the star. The magnitude of this tiny difference can 
be seen in the inset where the Kurucz model has been 
plotted with a dashed line. The photometric point at 0.36 
/im lies above the models because the Walraven B filter 
includes fiuxes on both sides of the Balmer jump, hence 
the average is higher than the local flux at the blue side 
of the jump. The agreement between the dereddened and 
synthetic photometry is compatible with the presence of a 
weak accretion shock signature in the ultraviolet colours 
of the star. 



|S1 and match the dereddened data (filled circles) so well 
that they are barely visible in the graph. In the inset we 
display an enlargement, in linear scales, of the region to 
the blue side of the Balmer jump to show the subtle dif- 
ference between the stellar (dashed line) and stellar plus 
shock (solid line) fluxes. 

We used the x^-minimization method to check if some 
UV excess was present in HD 34282. For three typical val- 
ues of the column energy fluxes (namely 1 x 10^"'^, 3 x 10^^ 
and 1 X 10^^ erg cm~^ s~^) we added the accretion flux 
to the Kurucz model with different filling factors / and 
performed the estimation for the fittings. Since the 
Walraven photometry was not contemporaneous with the 
EXPORT observations, the best-fit value for E{B—V) was 
0.08 in this case. The quality of the fits to the data is very 
good, but, as the observational errors are much larger than 
the difference between models with and without accretion 
flux, we can not distinguish between these models. This 
happened for all three values of the column energy flux 
used and in all cases the mass accretion rate derived was 
very similar, the average value being M= (8.2±0.2) x 10^^ 
M0 yr~^. This value is consistent with the low accretion 
rates expected in HAeBe stars (Hartmann et al. 1993; 
Ghandour et al. 1994; Natta et al. 2001) and also with 
estimates for UX Ori-like stars (Tambovtseva & Grinin 
2000). 

Given the tiny difference introduced by the shock and 
the small values obtained, it is clear that the excess flux 
is virtually impossible to measure in this case and that this 
determination constrains the mass accretion rate mostly 
as an upper limit. Once this method has been applied to 
other HAeBe stars with different amounts of IR excess, 
we will have better knowledge of its reliability. However, 
here we will adopt this value for the mass accretion rate 
toward the star because, being compatible with the ultra- 
violet photometry, it is also consistent with the physical 
parameters needed to fit the SED millimetre fluxes (see 
section 16.2.2(1 . In the next three subsections we describe 
in detail how three different components account for the 
observed SED of HD 34282. 



main-sequence value in cases with high mass accretion rate 
and large filling factors. 

Shock models for a central star with the HD 34282 
stellar parameters were computed. The energy flux of the 
accretion column, F, is a free parameter that controls the 
SED emitted by the accretion shock and the total mass 
accretion rate is proportional to this energy flux and to 
the filling factor, /, for the accretion column (MuzeroUe et 
al. 2003a). 

Fig. |S1 shows the fit of the photospheric plus accre- 
tion flux (solid line) to the dereddened Walraven photom- 
etry. Synthetic Walraven fluxes have been estimated by 
convolving the total flux (Kurucz plus shock fluxes) with 
the Walraven passbands, and compared to the dereddened 
data. These synthetic points are plotted as squares in Fig. 



6.2.2. The long wavelength part of the SED 

The flux at 1.3 mm is proportional to the mass of the disk 
since its thermal emission is nearly optically thin at this 
wavelength (see e.g.. Mannings 1994, Natta et al. 1997). 
Therefore we can use the flux density at 1.3 mm (Pietu et 
al. 2003) to estimate the viscosity parameter a of the sys- 
tem. Under the assumption of steady accretion, the disk 
surface density scales as E oc Ma~'^T~^ , S being the sur- 
face density, M the mass accretion rate, a the viscosity 
parameter, and T the temperature of the mid-plane disk 
(D'Alessio et al. 1999). For moderate accretion rates, the 
heating of the disk is dominated by stellar irradiation and 
T will not change very much with the mass accretion rate, 
so the mass of the disk will be roughly constant for con- 
stant Ma"^. Usual values of a range from 0.1 to 0.001, 
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0.01 being the one that reproduces best the milhmetre 
fluxes in T Tauri stars (D'Alessio et al 2001). We find 
that OL = 0.0008 yields the appropriate milhmetre flux for 
HD 34282, given the stellar and disk parameters consid- 
ered. 

The slope of the SED at millimetre wavelengths amm 
is normally modelled with the p dependence of the opac- 
ity Kx = kq{\/\o)~'^ (Beckwith et al. 1990; Beckwith & 
Sargent 1991) where amm — /? -I- 2 to account for the 
growth of the dust grains. In our case, the millimetre 
slope is controlled by the dust distribution and the op- 
tical depth of the disk. To obtain the observed millimetre 
slope anim = 2.2 (Fj, oc i^^-^), a combination of p = 2.5 and 
a maximum grain size flmax = 1 cm was needed. This value 
of p can be interpreted as the resulting size distribution of 
a coagulation process in which most of the disk dust mass 
resides in large grains (Miyake & Nakagawa 1993). 

The contribution of the mid-plane disk to the SED is 
shown in Fig. as a long-dashed line. There is an excess 
flux with respect to the observed flux at 3.4 mm measured 
by Pietu et al. (2003). We consider the sub-millimetre and 
millimetre slope well fitted with F^, oc v^'^ since five flux 
densities, between 0.45 and 2.6 mm, are consistent with 
it, and no significant change in the slope is expected (see, 
D'Alessio et al. 1990, but also Qi et al. 2003 for a counter- 
example) . Missing flux, a sudden change in the dust opac- 
ity at large wavelengths or the contribution of optically 
thick regions at small radii may explain the discrepancy 
(e.g. D'Alessio et al. 2003). A way to check the validity 
of this result would be to simulate the observation of this 
model for the mid-plane disk in the same line observed 
by Pietu et al. (2003), namely the ^^CO J = 2 ^ 1 tran- 
sition, and compare the resulting map with the observa- 
tions. This work is being done and it will be presented 
elsewhere. 

6.2.3. The mid-infrared region of the SED 

As can be seen in Fig. the IR emission coming from a 
single irradiated accretion disk model does not reproduce 
the observed mid-IR fluxes. According to the CIRR val- 
ues in the IRAS PSC, the cirrus infrared emission around 
this star is not an important contribution to the measured 
fluxes at 60 and 100 /im, so we consider that the flux in 
this part of the SED is also coming from the disk, implying 
that a second component is required. 

The IR excess produced by an irradiated accretion 
disk model with dust and gas well mixed decreases as 
the size of the dust particles increases (D'Alessio et al. 
2001). For a fixed dust mass, a decrease in the amount of 
large dust particles obviously implies the presence of more 
small grains in the disk; the effect is that the opacity of 
the dust to the stellar radiation increases and the height 
where most of the stellar radiation is deposited is greater. 
In that case, the disk intercepts more radiation from the 
star and hence it will emit more IR flux. In order to fit 
the mid-IR excess in HD 34282, we need the contribution 



of a disk component that we will call the surface disk with 
Omax = 1 Min, p ~ 3.5, the same mass accretion rate as 
the mid-plane disk (in this way we assure that the mass 
transfer is conserved through the disk) and a = 0.008. 
The contribution of the surface disk to the SED is shown 
in Fig. El as a dotted line. 

Considering both components, the disk would contain 
large grains in the mid-plane, which are responsible for 
the almost optically thin emission at sub-millimetre and 
millimetre wavelengths and a population of small grains - 
the surface disk - on top of the previous model, which is 
responsible for the emission in the mid-IR. The total mass 
of the disk around HD 34282 can now be estimated from 
the masses of both disk components (Table EJ assuming 
that the mid-plane disk extends out to some intermediate 
height and the surface disk occupies the higher layer to 
the end of the disk. Given that the disks are in hydro- 
static equilibrium, the density distribution with height is 
a gaussian function centered at the mid-plane. The total 
mass is then dominated by the mid-plane disk and lies 
between 0.6 and 0.7 Mq. The disk mass is less than half 
the stellar mass, supporting the hypothesis that the main 
gravitational field is due to the central star. 

Natta et al. (2000) used the flux at 2.6 millimetres to 
estimate a mass of 0.05 Mq for the disk of HD 34282, 
an order of magnitude smaller than our value. Differences 
in the gas temperature and the opacity parametrization 
used in their calculation may account for the difference. 
Nevertheless, the upper limit of the disk masses computed 
by Natta and collaborators for a large sample of HAeBe 
stars matches our estimate for the mass of the disk around 
HD 34282, showing that it is quite massive compared to 
typical HAeBe disks. 

6.2.4. The near-IR wavelength region of the disk 

Some HAeBe stars show conspicuous near-IR bumps. 
Hillenbrand et al. (1992) used high mass accretion rate 
models to fit their SEDs but Hartmann et al. (1993) ar- 
gued that such large values for M were incompatible with 
the bolometric luminosity of the objects. Finally Natta 
et al. (2001) proposed the presence of a frontally illumi- 
nated wall at the dust destruction radius where this near- 
IR bump would be emitted making compatible the shape 
of the SED and the luminosity of the object. 

We assume that the excess emission at 3 /im in the 
SED of HD 34282 originates in a wall-like structure at the 
inner part of the disk and model it as a black body at 
1400 K. Following Natta et al. (2001) and Dullemond et 
al. (2001), if this wall of irradiated material were optically 
thick then it would radiate as a black body at the local 
temperature, namely that of dust destruction (according 
to Pollack et al. 1994, rcvap~1200 - 1400 K). 

The position, i?rim, and height, iJrim, of this irradiated 
wall can be computed following the methods described in 
Dullemond et al. (2001). Also, MuzeroUe et al. (2003b) 
have recently fitted the near-IR bump in the SEDs of 
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Fig. 9. Disk models fitting the full SED of HD 34282. 
Filled triangles are dereddcned Walraven photometry 
from de Geus et al. (1990). Filled circles are dereddened 
EXPORT optical and near-IR photometry. Filled squares 
are LM magnitudes from Sylvester et al. (1996). Open 
symbols are the observed fluxes. The crosses are the IRAS 
colour-corrected fluxes and the open pentagons are sub- 
millimetre fluxes from Sylvester et al. (1996), Pietu et al. 
(2003) and Mannings & Sargent (2000). The dashed hne 
is the Kurucz model for the stellar parameters computed 
in Sectional The thin solid line is the 1400 K black body 
fit to the near-lR excess. The dotted line is the contribu- 
tion of the surface model with inclination angle i = 56°, 
M 8.2 X 10-9 Moyr-i, a = 0.008, i?i„==0.31 AU, 
-Rout=705 AU, p = 3.5 and Omax = 1 /ini. The long-dashed 
line is the contribution of the mid-plane disk model with 
the same parameters as the surface one except p = 2.5, 
Omax = 1 cm and a = 0.0008. The thick solid line on the 
top of the others is the addition of all components. See 
the text for explanation. 



Classical T Tauri stars with a black body arising from 
a similar model for a dust destruction wall. These authors 
find best-fit temperatures for the wall to be around 1400 
K. The main difference with respect to the first model 
is that the dust is not treated as optically thick but the 
position of the rim is computed with the same dust size 
distributions as in D'Alessio et al. (2001) and the methods 
used in Calvet et al. (1991, 1992). 

The projected area. A, of the rim is computed as in 
appendix B in DuUemond et al. (2001). The solid angle 
subtended by this area is A/cP, where d is the distance 
to the source. The normalization constant needed to fit a 
black body to the near-IR photometry is proportional to 
the solid angle subtended by the rim, so after computing 



the position of the wall, i?rim, we can estimate the cor- 
responding height so that the emitting area is consistent 
with the photometry. 

Then, using equation (14) in DuUemond et al. (2001) 
in a simple iterative procedure, we get i?rim = 0.31 AU 
and i?rim = 0.06 AU. Also, using the equation for 
in MuzeroUe et al. (2003b) similar results were obtained 
(namely i?rim = 0.31 AU and i7rim = 0.04 AU). This agree- 
ment is due to the fact that the different optical depth 
assumed for the dust in MuzeroUe et al. (2003) is com- 
pensated by the autoirradiation correction introduced by 
DuUemond et al. (2001). 

It is remarkable that the height obtained in this way 
is similar to the height of the surface disk component at 
^rim (-ff model = 0.06 AU) . However, we must take into 
account that the height of a single disk where the large 
grains are in the mid-plane and the small grains in upper 
layers may be different to this^. Firmer conclusions about 
these heights require a self-consistent determination of the 
vertical structure with settled dust, which is left for future 
work. 

As quoted in Natta et al. (2001), the hypothesis of 
the irradiated wall assumes a negligible opacity of the gas 
disk inside i?rim so that the stellar radiation is able to 
reach the inner wall. The calculation of the actual gas 
opacity inside the wall is difficult because of the multiple 
atomic and molecular absorption lines present. MuzeroUe 
et al. (2003b) studied the opacity of a gaseous disk inside 
the dust destruction radius using an extended set of gas 
opacities and a detailed treatment of the heating of the 
gas due to accretion and stellar irradiation. They found 
that for typical parameters of HAeBe stars, the optical 
depth of the inner hole would be low enough to allow a 
substantial fraction of the stellar radiation to reach the 
dust sublimation wall. 

The thin solid line in Fig. El is the contribution of the 
irradiated wall to the SED. The thick solid line is the 
addition of the Kurucz model and the contributions to 
the IR excess of the wall and the mid-plane and surface 
disks. In Table [S] a summary of the input parameters for 
the models is given. Also the values found for the height 
at 10 AU and the disks' masses for each component are 
given. 

6.3. The disk model for HD 141569 

The SED of this star has already been fitted with opti- 
cally thin dust models (e.g. Sylvester et al. 1996, Malfait 
et al. 1998). Recently Li & Lunine (2003) have published 
a highly detailed, thorough model in which the star is 
taken to be a Vega-type object with a debris disk filled 

This is a consequence of the fact that the model with large 
cimax has a smaller vertical scale height - is less flared - than 
the model with small flmax (see Table |3 , so the dust in the 
composite model would look like a slighly flared distribution of 
large particles surrounded by a more flared structure of smaller 
particles. 
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Table 5. Disk model parameters and results. 





HD 34282 


HD 141569 




Mid-plane 


Surface 




M (Mr vr"^) 


8 X 10"'-' 


8 X 10"'^ 


1 X 10"" 


a 


8 X 10"* 


8 X 10"^ 


1 X 10"^ 


Rin (AU) 


0.31 


0.31 


0.24 


Rout (AU) 


705 


705 


428 


P 


2.5 


3.5 


2.5 




1 cm 


1 fim 


1 mm 


Inclination 


56° 


56° 


51° 


Hio AU (AU) 


3.06 


4.13 


2.33 


Mdisk (Mq) 


0.70 


0.06 


6.3 X 10-"* 



with optically thin porous dust. We have shown in previ- 
ous sections, however, that HD 141569 is a PMS HAeBe 
star and not a Vega-type MS star, surrounded by a CS 
disk. This disk is clearly more evolved than usual CS disks 
detected around other HAeBe stars, but it cannot be con- 
sidered as a 'typical' debris disk, as indicated by the ra- 
tio between the IR excess and the bolometric luminosity 
and the non-negligible gas content and activity. Thus, the 
HD 141569 disk appears to be in an intermediate stage 
between 'first generation' PMS disks and 'second genera- 
tion' debris disks. We take this as a justification to carry 
out an exploratory study in which we fit the SED of HD 
141569 with PMS disk models and then compare the re- 
sults with those from the debris disks models, looking for 
convergent aspects. Note also that A-type Vega-like stars 
with debris disks are older systems than HD 141569 (Song 
et al. 2001). 

The observational results strongly indicate that there 
is no significant mass accretion. Thus, we assume that the 
disk is heated only by irradiation from the central star 
and that no viscous dissipation from the mass accretion 
onto the star is present, which translates in the models 
we are using into an extremely low mass accretion rate. 
The central star has the following parameters (SectionlS)): 
Teff = 10000 K, R^ = 1.70 Ro, M, = 2.00 Mq, = 25.77 
Lq and d= 108 pc, where the stellar radius was calculated 
from the stellar mass and the surface gravity. For the disk, 
we will use an inclination angle i = 51° (Weinberger et al. 
1999) and an outer radius i?out = 428 AU, i.e. the value 
of 400 AU given by Weinberger et al. (1999) scaled to our 
estimated distance of 108 pc. A change in the outer radius 
of no more than 50 - 60 AU will not affect substantially 
the final results. The inner radius is fitted to be i?in=0.24 
AU to account for the lack of near-IR excess in the SED. 

Fig. 1101 shows the results of our irradiated disk model 
compared to the SED of HD 141569. The peak at 7.7 fim 
is a PAH feature (see Li & Lunine 2003) that has not been 
modelled because it is not included in the dust model used 
(D'Alessio et al. 2001, Pollack et al. 1994). We find that 
a disk with a negligible mass accretion rate, M ^ 10~^^ 
Mq yr"-'^, and a = 0.001 (which translates into a dust 
and gas distribution of S oc i?"^ ° with i;(10AU) = 3 g 
cm"^), and a dust distribution with amax = 1 nim and 
p = 2.5 fits reasonably well the multiwavelength SED of 



SEDs of HD 34282 and HD 141569 



-8 



o 

-12 



-14 



' A IB. 

-1 A ^ 


- 


^ \ 




ill \ 




1,1 \ 




lll ^ 

lll 

III 

tF 




l|l 


\ \ 


j 

[ 


\ \ 
\ \ - 

\ \ 

\ \ ' 

\ \ 
S \ 
\ \ 
\ \ 
\ \ 
\ \ 
\ \ T 


HD 141569 


\ \ 

\ \ 



0.1 1 10 100 1000 



A(^m) 

Fig. 10. Disk model fitting the SED of HD 141569. Filled 
triangles are dereddened lUE data. Filled circles are 
dereddened EXPORT optical magnitudes. Filled hexagons 
are dereddened Ji/JT magnitudes from Malfait et al. (1998) 
and filled squares are LM magnitudes from Sylvester et 
al. (1996). The corresponding open symbols are the ob- 
served fluxes. The small open circles are the ISO-PHOT 
measurements. The crosses are the IRAS colour-corrected 
fluxes, the open pentagon is the submillimetre flux from 
Sylvester et al. (2001) and the arrow is the upper limit flux 
from Sylvester et al. (1996). The dashed line is the Kurucz 
model for the stellar parameters computed in Section \S\ 
The thick solid line is the disk model with inclination an- 
gle i = 5r, M = 10-" Mgyr-i, a = 0.001, i?in=0.24 AU, 
-Rout=428 AU, p = 2.5 and flmax = 1 nim. See text for 
details. 

HD 141569. The fractional IR luminosity from our model 
is 0.011, consistent with the measured value (Section O. 
The total mass (gas and dust) of our disk is 6.4 xlO~* 
Mq which is equivalent to 213 M0. The mass of the disk 
derived with this model is higher than the typical masses 
encountered in Vega- like disks (Thi et al. 2001) but clearly 
smaller than those found in HAeBe stars (Natta et al. 
1997), as expected given the evolutionary status of HD 
141569. TableElgives a summary of the input parameters 
for the model as well as the mass of the disk and the height 
at 10 AU^. 

^ Note that, whereas the disk of HD 141569 seems to be more 
evolved than that of HD 34282, the estimated age of HD 141569 
is less than that of HD 34282 (see Table |^. However, one has 
to be careful when comparing this parameter for both stars, 
given the fact that HD 141569 is more massive and has larger 
metal abundance than HD 34282. Therefore, the evolutionary 
time scales are different for both objects. The initial conditions 
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In the model by Li & Lunine (2003) the surface den- 
sity of dust particles is introduced with three analytical 
approximations that reproduce the surface brightness pro- 
files observed in thermal and scattered light images by 
Augereau et al. (1999), Weinberger et al. (1999), Mouillet 
et al. (2001) and Marsh et al. (2002). However, the verti- 
cal height of the dust is computed assuming vertical hy- 
drostatic equilibrium with a gas temperature taken to be 
a power law of the distance from the central star. In our 
disk models, both the density profile and the vertical scale 
height are computed taking into account the stellar irradi- 
ation heating (assuming that dust and gas are in thermal 
equilibrium). Concerning the continuum emission in the 
infrared to millimetre wavelengths, the presence of gaps of 
moderate width with smaller dust density does not change 
dramatically the shape of the SED (see e.g. Fig. 3 in Li & 
Lunine (2003) where the three dust components are added 
to yield the smooth observed SED). Therefore we can as- 
sume that our continuous dust distribution is a good first 
approximation to the actual dust density with the advan- 
tage that it is self-consistently computed, along with its 
vertical hydrostatic structure, with the set of stellar pa- 
rameters previously determined. This is in agreement with 
the results obtained by Clampin et al. (2003). 

Concerning the dust, Li & Lunine (2003) found that 
a 'cold-coagulation' dust composed of porous dust grains 
with a porosity P — 0.9 (90 % of the volume is vacuum 
and 10 % material), and a power law distribution of sizes 
n{a) oc a~^'^ with 1 ^m < a < I cm, fitted very accu- 
rately the whole continuum SED of HD 141569. This dust 
model is not very far from our assumed n{a) oc a~^'^ with 
0.005 fim < a < 1 mm mostly at large sizes. The differ- 
ence in the optical properties of the compact dust grains 
compared to the porous dust grains may be accounted for 
by using a different power law exponent for the size distri- 
bution. The contribution of the PAH molecules to the IR 
emission is more important around the peaks at 3.3, 6.6, 
7.7, 11.3 ^m (see Li & Draine 2001 for details) and the 
continuum is well reproduced by the porous or compact 
dust populations. 

The main effect of having a disk with gas and dust, in- 
stead of a pure dusty disk, is that the heating due to stel- 
lar irradiation produces an increase in the vertical scale 
height of the disk, yielding a more dilute dust population. 
Summarizing, we find that our approach fits the observed 
SED of HD 141569 reasonably well, although the model 
we have used does not include as many details as that 
developed by Li & Lunine, particulary in the treatment 
of the dust properties. In this sense, our result is inter- 
esting since it shows the possibility of studying CS disks 
in transition objects from both perspectives, namely the 
PMS and MS sides. 



of the clouds where the stars and disks formed must also play 
a role in the evolutionary scenario. 



7. Conclusions 

We have presented a study of the properties and spectral 
energy distributions of the HAeBe stars HD 34282 and 
HD 141569 and their disks, based mainly on observations 
made by the EXPORT consortium (Eiroa et al. 2000). The 
main conclusions can be summarized as follows: 

- The spectroscopic and photometric observations of 
both stars have been described and analysed in detail 
with the aim of throwing light on their evolutionary 
status and the dynamics and variability of their pho- 
tospheres and disks. 

- Both stars are metal-deficient. Our estimates are 
[Fe/H]=-0.8 for HD 34282 and -0.5 for HD 141569. 

- The analysis has provided values for the stellar dis- 
tances. The distance to HD 34282, namely 348 pc, is 
consistent with the value given by Pietu et al. (2003) 
and corrects the Hipparcos distance that is uncertain 
due to the large error in the parallax. The new dis- 
tance for HD 34282, together with its low metallic- 
ity, resolves the anomalous position of this star in the 
HR diagram previously reported under the assump- 
tion that the star has solar abundances. The distance 
to HD 141569, namely 108 pc, is consistent with the 
Hipparcos data and in this case, the question of the 
anomalous position of this star in the HR diagram is 
resolved by using the correct set of tracks with the 
appropriate metallicity. 

- In addition to the above two parameters, our analy- 
sis has provided for both stars values for the effective 
temperatures, spectral types, masses, gravities, lumi- 
nosities, ages, projected rotational velocities and ex- 
tinctions. These results are given in Table 21 

- Complete SEDs from the ultraviolet to the millimetre 
range have been constructed for both stars and their 
disks using the EXPORT data plus additional results 
gathered from the literature and on-line databases. 
Table |31 gives a full account of this compilation and 
the data sources. 

- The self-consistent irradiated accretion disk models of 
D'Alessio et al. (1998, 1999 and 2001) have been used 
to fit the SEDs of both systems. A summary of the 
disk parameters obtained is given in Table [S] 

- The SED of the disk of HD 34282 has been repro- 
duced using a three-component model. The mid-plane 
disk is modelled with a maximum dust grain size of 1 
cm and emits mainly at submillimetre and millimetre 
wavelengths, the surface disk model has a maximum 
grain size of 1 fim and is responsible for the mid-IR 
excess emission. A near-IR bump similar to those seen 
in other HAeBe stars (see Natta et al. 2001) is fit- 
ted with the emission of a frontally illuminated wall 
placed at the dust destruction radius (DuUemond et 
al. 2001; MuzcroUe et al. 2003b, D'Alessio et al. 2003). 
The geometry of the wall is constrained by the near-IR 
photometry. 

- In the process of fitting the SED of HD 34282, we have 
tried for the first time a new method to estimate the 
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mass accretion rate towards a HAeBe star using ultra- 
violet Walraven photometry and the magnetospheric 
accretion shock models of MuzcroUc ct al. (2003a). Wc 
obtained un upper limit for the mass accretion rate 
in HD 34282 compatible with the parameters needed 
to fit the millimetre SED. The low value found does 
not support the idea by Herbst & Shevchenko (1999) 
about the FU Ori-like mechanism causing the photo- 
metric activity of UX Ori stars. 
- The SED of HD 141569 has been fitted with a disk 
model irradiated by the central star with negligible 
accretion and with large grains (a,„ax = 1 mm). A 
slightly flared dust and gas distribution was obtained, 
which is consistent with a substantial amount of resid- 
ual gas in the disk. 
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